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Abstract 

We discuss exclusive vector meson electroproduction within the QCD 
collinear factorization framework. In Bjorken kinematics the ampli- 
tude factorizes in a convolution of the nonperturbative meson distri- 
bution amplitude and the generalized partem densities with the per- 
turbatively calculable hard-scattering amplitudes, which are presently 
known to next-to-leading order (NLO). At small xb NLO corrections 
are very large. It is related to appearance of BFKL type logarithms in 
the hard-scattering amplitudes, that calls for a resummation of these 
effects at higher orders. Here we report the first results of such resum- 
mation. 

1 Introcuction 

The process of elastic vector meson electroproduction on a nucleon, 

l*(q)N(p)^V(q')N(p>), (1) 

where V = p°, u, <f>, was studied in many fix target and in HERA collider experiments. On 
the theoretical side, the large negative virtuality of the photon, q 2 = —Q 2 , provides a hard 
scale for the process which justifies the application of QCD factorization methods that allow 
to separate the contributions to the amplitude coming from different scales. The factorization 
theorem [1] states that in a scaling limit, Q 2 — > oo and xb = Q 2 /2(p ■ q) fixed, a vector meson 
is produced in the longitudinally polarized state by the longitudinally polarized photon and that 
the amplitude of the process (Q]) is given by a convolution of the perturbatively calculable hard- 
scattering amplitudes C\ the nonperturbative meson distribution amplitude (DA) <j)y, and the 
generalized parton densities (GPDs) IP. 

A=Y^ J J dx dz H l (x,£,t,(j, F ) C\x,z,fj, F ) (pv(z,fi F ) , (2) 
i=g,g 

where £ = xb/(2 — %b) is the skewness variable, t = (p — p') 2 and /j, F is a factorization 
scale. GPDs encode important information on hadron structure, including aspects that cannot 
be deduced directly from experiment, like the transverse spatial distribution of partons and their 
orbital angular momentum, for more details see [2]. 

Deeply virtual Compton process (DVCS) provides the theoretically cleanest access to 
GPDs. Recently two-loop effects were incorporated into the analysis of DVCS [3]. A theo- 
retical description of exclusive meson production is more involved since it includes an additional 



nonperturbative quantity, a meson DA. The primary motivation for the strong interest in this pro- 
cess (and in the similar process of heavy quarkonium production) is that it can serve to constrain 
the gluon density in a nucleon. Indeed, in vector meson production case the gluon GPD enters 
the description already at the leading order (LO) in the strong coupling a s , whereas in DVCS it 
appears first at NLO, and, like in inclusive DIS, is accessible only through scaling violation. 



2 NLO corrections 

The hard-scattering amplitudes for process (fl]) were calculated at NLO in [4], and for exclusive 
heavy quarkonium photoproduction in [5]. The analysis of NLO effects showed that in kine- 
matics typical for the HERA collider experiments, 15 ~ 1CT 3 , the NLO corrections are huge 
even for really large values of hard scales ~ 30 GeV 2 . If the factorization scale is chosen close 
to the value of a hard scale, /ip ~ Q, the corrections have opposite signs in comparison to 
the Born term. Which may lead to the change of signs of the imaginary and the real parts of 
the amplitude within phenomenologically relevant interval of xb- Besides, the factorization and 
renormalization scale uncertainties were found being very large. 

Recently these findings were confirmed in [6], where very detailed analysis of the cross 
sections and the transverse target polarization asymmetries in exclusive meson production was 
performed both for small and larger values of xb, typical for fixed-target experiments. For the 
fixed target kinematics it seems that NLO corrections start to be under control, though their 
values are still large at presently available values of Q 2 . For the transverse target polarization 
asymmetries the situation is better, in some cases. 

Going back to small xb, why NLO corrections are large in this case? The inspection of 
NLO hard-scattering amplitudes shows that the imaginary part of the amplitude dominates and 
that the leading contribution to the NLO correction originates from the broad integration region 
£ <C x <C 1, where the gluonic part approximates (N c = 3 is a number of colors) 
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Given the behavior of the gluon GPD at small x, H 9 (x, £) ~ xg(x) ~ const, we see that NLO 
correction is parametrically large, ~ ln(l/£), and negative unless one chooses the value of the 
factorization scale sufficiently lower than the kinematic scale. For the asymptotic form of meson 
DA, 4>y(z) = 6z(l — z), the last term in ([3]> changes the sign at fip = Q, for the DA with a more 
broad shape this happens at even lower values of up. Similar, m(l/£) enhanced, contribution 
appears also in the quark singlet channel. 

The partonic momentum fraction x is related to the Mandelstam energy variable s of the 
partonic subprocess x/£ ~ s/Q 2 . The leading part of NLO partonic amplitude (proper normal- 
ized) grows as the first power of energy, x ~ s, whereas at LO partonic amplitude behaves like 
a constant at large s. The reason for this difference is the appearance, stalling from NLO, of 
partonic diagrams with the gluon exchange in the t— channel, see Fig. 1. At LO one has only 
diagrams with the quark exchange, both for the gluon and quark channels. 

At higher orders the diagrams with gluon t— channel exchange give contributions to the 
amplitudes of partonic subprocesses enhanced, for n loops, as a" log" -1 x/x. In its turn, these 




Fig. 1: NLO diagrams with t— channel gluon exchange; the gluon and the quark GPD contributions. 

terms inserted in the factorization formula will produce large contributions ~ a™ log n (l/£) to 
the process amplitude, where each power of the strong coupling is compensated by the same 
power of a large logarithm of energy. It is a natural idea to resum these enhanced at small x b 
contributions using the BFKL approach [7]. 



3 High energy resummation 

The central point in this high energy resummation is to perform it consistently, without spoiling 
the all-order factorization of collinear singularities. A care should be taken of the factoriza- 
tion scheme used at the factorization of the process amplitude (O in terms of GPDs and hard- 
scattering amplitudes. The higher order terms of the hard-scattering amplitudes derived within 
the high energy approximation (BFKL approach) can be supplemented by the knowledge of 
hard-scattering amplitudes calculated exactly at fixed order. Then, one can use them together in 
factorization formula (O without double counting. 

For inclusive hard processes, heavy quark production and DIS, the method of such high 
energy resummation was elaborated in [8]. It is based on Curci, Furmanski and Petronzio ap- 
proach [9] to separation of collinear singularities. The amplitudes on a parton (quark, gluon) 
target are considered in D = 4 + 2e non-integer dimentions. That separates automatically the 
leading twist. Collinear singularities appear in this approach as 1/ e n poles, these poles are ab- 
sorbed into a definition of parton densities. Another essential ingredients of [8] method is an 
analysis of Mellin moments, high-energy terms in Mellin moment space N look like singulari- 
ties (a s /N) n at N — > 0, and a consideration of BFKL equation in D = 4 + 2e dimensions. 

We found3 that this technique may be directly generalized on the analysis of exclusive 
non-forward reactions. Below we present the first results of this study. 

Like in DIS the imaginary part of the a amplitude is given by the sum of quark singlet and 
gluon contributions 

l 

2mA(£, t) = - [ dx 

e 

£>W and are the imaginary parts of the quark and gluon hard-scattering amplitudes. In 
difference to forward DIS case the parton densities in (0]) depend on both longitudinal momentum 
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fractions. Due to that the Mellin moments of the amplitude do not factorize into the product of 
the moments 

D N (t) = J dU N lrnA(t,t) = (5) 
o 

1 1 r 

JJdudx u N - 1 x N £>(+) («) 7i {+) (x, ux,t) + ±D 9 (u) H 9 (x, u x, t) 

L 

Using polynomiality property of GPDs, in particular for the gluon case 



/ dxx n H 9 (x,rj,t) = £ (2riyA° n+2ij (t) + (2r,r +2 C 9 n+2 (t) 

J ■ n 



(6) 



j =0, even 

one can show that for the integer odd N 



D N (t) = ]T 2 fc [^5S*-i^4*+^(*) + D 9 N+k _ 2 A% +k+hk (t)] . (7) 

k=0 

Which is a sum of moment products (not just a product, as in DIS case). 

One can analytically continue (O from the integer odd N into entire complex N plane. The 
high energy asymptotic of the amplitude is related with the behavior of -Djv(i) near unphysical 
point N — > 0. One can split © into a sum of the singular and the regular at N — > parts 

D N (t) = cfrkftHt) + C 9 N g N (t) + D r ^(t) (8) 
The singularities of the sum f7]) at N — > are due to k = term only. Therefore 

4 +) =4T-i> ?£ ) (*) = ^W*). ^ = ^-2. 5iv(t) = ^ +li o(t). (9) 

Note that at t — > 0, q^\t) gw{t) reduce to the moments of usual parton densities 

l l 
?#°(t) 4 +) = / <feA (+) (i) , ffiv(t) - 97V = / ^/^W • (10) 



This consideration shows that a non-forward nature of hard exclusive reactions does not 
complicate much their analysis in the high energy limit. Therefore the method used in DIS [8] 
may be applied here. The difference between DIS and our case is in the different form of k t 
dependent amplitudes for corresponding partonic subprocesses. 

Below I will concentrate on the dominant at high energy gluon contribution. The results 
will be presented for the process (Q]) ( assuming for simplicity the asymptotic form of meson DA) 
and for the process of heavy quarkonium electroproduction (where the formation of quarkonium 
is treated in NRQCD). The amplitude is presented as follows 

l 

ImA°~H°(Z,Z)+ f -ff»(^)^C„(L)-^-lor^ (11) 
J x In — 1)1 £ 
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here a s = N c a s /Ti, we omitted normalization factors irrelevant for the subsequent discussion, in 
the r.h.s H 9 (^, £) represents the Born contribution and the sum stands for the high energy terms. 
C n (L) are the polynomials of variable L = log which we need to calculate. 

Note that the Born term belongs to the regular part (in terms of (0), whereas the high 
energy terms behave as (a s /N) n at N — > 0. Therefore in the high energy terms one can replace 
gluon GPD in (fTTT > by its forward limit, H 9 (x, £) — > xg{x), but in the Born contribution H 9 (£, £) 
should be kept different from xg(x). 

Omitting all details of the derivation we just present the results. We work in MS scheme. 
We define (properly normalized) k t dependent amplitude of the gluon subprocess 

M»?) - / * v + , )q , *M//* ^ , d2) 

then we calculate its Mellin transform 

oo 2 / o \ 

dk t I k t 




M7) = 7 / I I h v (kt) = L r[2 ^ 7] " ■ (13) 



o 

The high energy terms are defined from the expression 

C%~h v (i)R(^\ ■ (14) 

The gluon anomalous dimension is determined by the solution of equation 1 = {a s /N)x{^i), 
where x(l) i s trie BFKL eigenfunction, function R depends on a s /N and is defined in [8]. 
Expanding C 9 N in the series of variable y = a s /N one can obtain analytical expressions for the 
polynomials C n (L). 

Below we illustrate the values of these polynomials for the case fj,p = Q 

V : l-2y + 4y 2 -2.39y 3 -4.09y 4 + ... 
onium : 1 - 1.39 y + 2.61 y 2 + 0.481 y 3 - 4.96 y 4 + . . . 
F L : 1 - 0.33y + 2.13y 2 + 2.27 y 3 + 0.434 y 4 + ... 

here the first two lines represent results for the exclusive light vector meson and quarkonium 
production respectively, in the third line we show for comparison the results for longitudinal DIS 
structure function [8]. We see that numerical values of Ci(0) are negative in all case, but for 
the exclusive reactions its absolute values are about 4-^6 times larger then in the case of Fl, 
explaining very large negative NLO corrections found for exclusive meson production. On the 
other hand, the values of the second polynomial are positive and large, 62(0) = 4 for the light 
vector meson production. This gives a hope that inclusion of these high energy terms in the 
analysis may stabilize predictions for exclusive meson production. 

To investigate this possibility we perform the following numerical study. We calculate the 
amplitude of light vector meson production with (fTTT) . where in the high energy terms we use 



very simple input for the gluon density H 9 (x, £) ~ xg(x) ~ x~ 0,2 , for the Born term we take 
H 9 {^ : £) = 1.2 £ <?(£). Definitely, more realistic input for gluon GPD should be used (especially 
for H 9 (£,£)), but at the present stage we just want to clarify the qualitative role of the high 
energy terms. In Fig. 2 we present the energy dependence of the amplitude (in arbitrary units) 
calculated for two values of photon virtuality Q 2 = 10 and 20 GeV 2 , for the running coupling 
we use a s (10) = 0.25 and a s (20) = 0.16, and for the factorization scale fj,p = Q 2 /2. The solid 
line on Fig. 2 represents the Born contribution, the dashed line - the Born + the first high energy 
term, the dotted line - the Born + 2 first high energy terms, the dashed-dotted - the Born + 6 
first high energy terms. We see that high energy resummation is convergent fast, the difference 




Fig. 2: The convergence of the high energy resummation: (a) Q 2 = lOGeV 2 , (b) Q 2 = 20GeV 2 . 

between the dashed and the dashed-dotted lines is not big for Q 2 = 10 GeV 2 and is really small 
for 20 GeV 2 cases. The other observation is that the inclusion of only first high energy term 
(dashed line) seems to be a bad approximation. Even for 20 GeV 2 case, where the Born and the 
resummed results are close to each other, the dashed line is about factor of 3 below. 

The dependence of the amplitude on the choice of factorization scale is shown in Fig. 3. 
Again, the dashed lines correspond to the Born + the first high energy term, the dashed-dotted 
lines - the Born + 6 first high energy terms. The upper dashed and dashed-dotted lines are 
for n F = (J 2 /4, the lower dashed and dashed-dotted lines are for /j, F = Q 2 /2. We observe 
sizable reduction of the factorization scale dependence if the high energy terms are resummed in 
comparison to the case when only the first of these terms is taken into account. 

4 Summary 

Large NLO corrections are found for hard exclusive vector meson production. At intermediate 
to larger values of xb, typical for fixed-target experiments, it seems that NLO corrections start 
to be under control for the large values of Q 2 , say above 10 GeV 2 . However, the situation 
is much worse for the region of small xb, typical for the HERA collider experiments. Here 
NLO corrections are not under control even for such large values of hard scales as 30 GeV 2 , 
which prevents the interpretation of the precise HERA data in terms of GPDs. The problem is 
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Fig. 3: (a) Q 2 = lOGeV 2 , (b) Q 2 = 20GeV 2 



related to appearance of BFKL type logarithms in the hard-scattering amplitudes, that calls for a 
resummation of these effects at higher orders. Here we present the first results for such study. The 
methods used earlier for forward DIS process may be generalized to the case of nonforward hard 
exclusive reactions. We obtained analitical results for the corresponding high energy terms in 
(fTTb . The first numerical calculation incorporating the high energy resummation is encouraging. 
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